DGMRES: A GMRES-type algorithm for Drazin-inverse solution of singular nonsymmetric linear systems
نویسنده
چکیده
In a recent work by the author [Linear Algebra Appl. 298 (1999) 99] Krylov subspace methods were derived for the Drazin-inverse solution of consistent or inconsistent linear systems of the form Ax = b, where A ∈ CN×N is a singular and in general non-hermitian matrix that has an arbitrary index. One of these methods, modeled after the generalized conjugate residual method (GCR) and denoted DGCR, is considered in the present work again. It is shown that all of the approximations produced by DGCR exist, and a GMRES-like algorithm, denoted DGMRES, for its implementation is derived. Like GMRES, DGMRES too is economical computationally and storagewise. © 2001 Elsevier Science Inc. All rights reserved. AMS classification: 15A06; 15A09; 65F10; 65F50
منابع مشابه
GGMRES: A GMRES--type algorithm for solving singular linear equations with index one
In this paper, an algorithm based on the Drazin generalized conjugate residual (DGMRES) algorithm is proposed for computing the group-inverse solution of singular linear equations with index one. Numerical experiments show that the resulting group-inverse solution is reasonably accurate and its computation time is significantly less than that of group-inverse solution obtained by the DGMRES alg...
متن کاملSingular constrained linear systems
In the linear system Ax = b the points x are sometimes constrained to lie in a given subspace S of column space of A. Drazin inverse for any singular or nonsingular matrix, exist and is unique. In this paper, the singular consistent or inconsistent constrained linear systems are introduced and the effect of Drazin inverse in solving such systems is investigated. Constrained linear system arise ...
متن کاملA Bi-CG type iterative method for Drazin-inverse solution of singular inconsistent nonsymmetric linear systems of arbitrary index
متن کامل
Optimal Control for Descriptor Systems: Tracking Problem (RESEARCH NOTE)
Singular systems have been studied extensively during the last two decades due Abstract to their many practical applications. Such systems possess numerous properties not shared by the well-known state variable systems. This paper considers the linear tracking problem for the continuous-time singular systems. The Hamilton-Jacobi theory is used in order to compute the optimal control and associa...
متن کاملAdditional results on index splittings for Drazin inverse solutions of singular linear systems
Given an n × n singular matrix A with Ind(A) = k, an index splitting of A is one of the form A = U − V , where R(U) = R(Ak) and N(U) = N(Ak). This splitting, introduced by the first author, generalizes the proper splitting proposed by Berman and Plemmons. Regarding singular systems Au = f , the first author has shown that the iterations u(i+1) = U#V u(i) + U#f converge to ADf , the Drazin inver...
متن کامل